Analytical results on the Muller's ratchet effect in growing populations 
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Fontanari et al introduced [Phys. Rev. Lett. 91, 218101 (2003)] a model for studying the 
Muller's ratchet phenomenon in growing asexual populations. They studied two situations, either 
including or not a death probability for each newborn, but were able to find analytical (recursive) 
expressions only in the no-decay case. In this paper a branching process formalism is used to find 
recorrence equations that generalize the analytical results of the original paper besides confirming 
the interesting effects their simulations revealed. 
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It is widely recognized that the rate of deleterious mu- 
tations being much higher than that of either reverse or 
beneficial mutations can be a serious threat to the sur- 
vival of populations at the molecular level. About forty 
years ago, H. J. Muller conjectured [l| that, in such con- 
ditions, the mean fitness of finite lineages lacking mech- 
anisms of genetic repair should decay with time, due to 
the sucessive loss of the fittest individuals. This stepwise 
fluctuation-induced phenomenon is known as Muller's 
ratchet [2| and have received growing attention in re- 
cent years. It has been argued to be responsible for the 
origin of some diseases, for the fitness decay in some ex- 
periments with microorganisms and even for the origin 
of sex, as summarized in [l, [j] , although many of these 
claims still lack conclusive evidence. 

In contrast with the first studies [f|, recent work on 
Muller's ratchet have focused on models with variable 
population size, either under growingconditions 0, 0] or 
mixing growth with bottlenecks [1, Q. Indeed, this is 
the realistic condition that holds both in nature and in 
controlled experiments [1 01 ] . with exponential growth of 
microorganisms and bottleneck transfers. 

Specifically, in 0] Fontanari and collaborators (hence- 
forth FCH) introduced a fully stochastic discrete time 
model for growing asexual lineages in which each mem- 
ber of the population is replaced in the next generation 
by a number k of descendants distributed according to a 
Poisson distribution of parameter A = R, i.e., with prob- 
ability 



p{k; A) = exp (—A) 



A^ 
fe! 



(1) 



Each newborn acquires a number of new mutations also 
Poisson distributed, but now with mean value U . In 
addition, an individual with j mutations has a chance 
to leave offspring only if it survives, what happens with 
probability (1 — s) 3 ', where s G [0, 1] is a selective coeffi- 
cient. FCH focused on finding conditions for the halting 
of Muller's ratchet when the population is founded by a 



single mutation- free individual (called master sequence). 
They used three recursion equations to study the neu- 
tral case s = and recurred to simulations in the case of 
evolution under decay. 

The aim of this Brief Report is to show that the the- 
ory of branching stochastic processes, both in its simple 
(just one type of individual) and multi-type versions, al- 
lows a complete description of FCH's model, not only 
recovering their analytical results in the neutral case in a 
straightforward way, but also giving expressions valid for 
the case s > 0. In particular, it is shown that the coun- 
terintuitive acceleration of the ratchet's activity with an 
increase on the selective coefficient found through simu- 
lations in Q is a real phenomenon and not an artifact. 
The theory of branching processes is a veritable subject, 
originally developed by J. B. S. Haldane [HI and R. A. 
Fisher [lj, [H| m the birth of the modern population ge- 



netics and today is described in many textbooks like [Fj 
The multi-type theory have found lots of applications in 
evolutionary dynamics, e.g. [l5l |. 

To start with, we consider the neutral case. Defining 
M{k,t) as the probability that the population is com- 
posed of k individuals in generation t, regardless their 
mutational load, FCH found that the generating func- 
tion g(z, t) = X^fcLo z N{k, t) obeys the equation 



g(z,t+ 1) = g[e 



-fl(l-z) 



(2) 



Being the probability that the population is not ex- 
tinct at generation t, = 1 — g(0,t) and they argued 
that 



S( t+1 ) = l-exp(-RS®). 



(3) 



In a general simple branching process (henceforth SBP), 
if there is a single founder and (p(z) is the generat- 
ing function of the progeny distribution f(k), tp(z) = 
££o *"/(*). then 



g(z, t+l) = g[(p(z), t] = ip[g(z, t)] 



(4) 
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and it is easily seen that Eqs. ([2]) and ([3]) are just par- 
ticular instances of the above equation in the Poissonian 
case, f(k) = p(k; R). 
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It is another standard result in the theory of the SBP 
that the mean number of offspring of an individual been 
strictly greater than 1 is a necessary and sufficient con- 
dition for the asymptotic survival probability of the pop- 
ulation to be greater than zero. In the neutral case, ne- 
glecting the genotype of all individuals is harmless: they 
can be considered all equivalent. Thus, the evolution of 
the system is properly described by a SBP and, when 
s = 0, a necessary condition for the survival of the pop- 
ulation is R > 1, a result FCH derived from Eq. ([3]). 

More results follow from a decomposition theorem [l6| 
that is essential to this work. Roughly speaking, it states 
that if each "object" generated by a Poisson process with 
mean R is independently attributed to a class with prob- 
ability pi, then the number of objects in class i is again 
Poisson distributed, with parameter Rpi, and is indepen- 
dent of other classes. This "allocation" is exactly what 
mutation does and the classes are defined by the muta- 
tional load. 

Being e~ u the probability of a perfect replication, each 
individual originates an average of Re~ u error-free repli- 
cas of itself, following a Poissonian SBP. Therefore, if 
Re~ u < 1, the mutation-free subpopulation certainly 
goes extinct. Remembering we still are in the no-decay 
case, this argument applies qualitatively to all classes of 
individuals, even whether a new class that had just be- 
came the least-loaded one has more than one member at 
that moment. Hence, if R > 1 but Re~ u < 1, all classes 
of sequences have a finite lifetime and Muller's ratchet 
never halts in the surviving populations. Conversely, if 
Re~ u > 1, there is a positive probability (S(°°) that the 
least-loaded class survives the expansion process, given 
by the asymptotic solution of Eq. ((3]) when R is replaced 
by Re~ u . Every time the ratchet clicks is equivalent to 
a new start of the process, always with a positive sur- 
vival probability. Hence, if the population does not get 
extinct, some class will survive indefinitely and halt the 
ratchet. 

FCH noticed the criterium they found for the halting 
of Muller's ratchet (U <h\R or, equivalently, Re~ u > 
1) implies that, on average, each master sequence must 
generate more than one perfect replica in order to the 
population to be viable. However, despite recognizing 
that criterium as intuitive, they could not predict it in 
advance. Actually, they derived it from a recursion for 
V n (t), the probability that at generation t the minimum 
number of mutations in the population be n, valid only 
in the neutral case. It turns out to be the most relevant 
quantity for the study of Muller's ratchet since it reveals 
the time dependence of the mean number of mutations of 
the fittest member of the population. So it is noteworthy 
that the theory of a multi-type branching process allowed 
the finding of the main result of this paper, an analogous 
recursion for V n (t) in the general case s 6 [0, 1], which 
will be discussed from now on. To be honest, there is an 
infinite number of types of individuals, since there is no 
such thing as a maximal mutational load. But the results 
inspire confidence in the adopted approach. 



At this point, some remarks on notation are necessary. 
A vector u has infinite components and the first index 
is 0, to account for the mutation-free class. Explicitly, 
a phrase like "the population is in state k" means that 
k = (ko, ki, . . .) and there are kj individuals in class j. 
It is also convenient to define u = (u' 3 ',uW), where 
u' 3 ' = (ito, . . . , Uj-i) and uW = (ibj, Uj+i, . . .) for any 
j > 1. The same rules apply to constant vectors too. So, 
0^) means a vector with just j components, all null, and, 
despite how close lw appears to be to 1, they are not 
the same object, since the first index referred to in lW) 
is j. 

Let /i(k) be the probability of an i-individual to gen- 
erate offspring k and Mi (k, t) be the probability that the 
population is in state k at generation t, given a single 
founder with i mutations. The generating functions as- 
sociated to these two joint distributions are 

oo 

^( Z )=^/ i (k)n#, (5) 

k j=0 

and 

oo 

^(z,t)=^M(k,t)n^'. ( 6 ) 

k 3=0 

respectively. The generalization of Eq. (j4]) to the multi- 
type setting is 

g(z,i)= V [g(z,<-l)], (7) 

when there is only one founder. Thus, to characterize 
the population at a given time, whatever be the initial 
condition, it is necessary to know its properties in the 
preceding generation as if it had originated from all pos- 
sible types of founders. 

Inquiries about extinction acquire a broder sense in 
this case, since now it is possible to talk about the sur- 
vival of specific classes inside the community. Let qi(j,t) 
be the probability that the smallest index of a populated 
class at time t be at least j, given a founder in class i. 
From Eqs. © and (0), 

qtf.f) = g[(0 (i) ,l«),t] = ¥>[q(j,t- 1)]. (8) 

It is clear that Vu it) , the probability distribution of the 
smallest index j of a class still alive at t, given a founder 
in class i, is given by 

V ij (t) = q i {j,t)-q i {j + l,t), (9) 

and, since V n it) — Vo n {t), this is all the information 
needed to study Muller's ratchet. 

The existence of analytical expressions for the multi- 
dimensional generating functions is essential for this pro- 
posal to be useful. At this point, the decomposition the- 
orem of Poisson processes enters again. It is important 
to notice that, under decay, each mutant is temporar- 
ily allocated on some class, depending on its mutational 



load, but after that it may not survive and thus may be 
"redirected" to a "sink class" that plays no role on the dy- 
namics. Therefore, the factor pi informally introduced in 
the presentation of the decomposition theorem must take 
into account both mutation and the survival probability, 
while the number of individuals in any class still is Pois- 
son distributed with parameter Rpi. For instance, the 
number of individuals with i mutations directly descend- 
ing from a j-mutant is given by a Poisson distribution 
with mean Rp(i — j; — s)\ Of course, independence 
of classes still holds and this allows a convenient fac- 
torization of the joint distribution /^(k) as the product 
of Poisson distributions. Consequently ^(k) also fac- 
torizes (as a product of Poissonian generating functions) 
and, after some algebra, it follows from Eq. ([8]) that 

j'-i 

uU, t) = n ex p i R p( k - *; u ) (! - s ) k ■ 

k—i 

■[q k (j,t -1)-1]}. (10) 

The product is finite because qi{j,t) = 1 when j < i. 
Some time after discovering this result, we were informed 
that it is a special case of a general theory developed in 
[HI and [j"?]]. Nonetheless, its present application to the 
study of Muller's ratchet in FCH's model is still unre- 
ported. 

When s = 0, Qi(j,t) depends only on j — i and Eqs. 
(JH) and (fTOf together are equivalent to the recursion for 
V n (t) derived in 0], as expected. And, under decay, they 
still can be iterated easily. As an example, Fig. [T] (similar 
to Fig. 2 in Q) illustrates an unexpected phenomenon 
discovered by FCH: at a given generation, it is possi- 
ble that the average mutation load of the least-loaded 
class increases if selection gets stronger, R and U kept 
fixed. Indeed, this effect seems counterintuitive at first. 
But since fitness is absolute in the model, any expla- 
nation must not rely on competition in the population 
and the analysis may not be straightforward. It seems 
valid for the author to reverse the reasoning: why should 
an increase in the intensity of decay always imply a de- 
crease in the average mutational load of the fittest class 
in such a complex model, in which fitness is absolute and 
any change in parameters affect the probability of extinc- 
tion? It is important to remember that all averages are 
calculated with probabilities conditionated on the sur- 
vival of the population. Maybe the enhanced activity of 
Muller's ratchet seen in Fig. Q] be just a consequence 
of some peculiar dynamics of the spectrum of the pop- 
ulation (regarding the distribution of individuals among 
the mutational classes) in some region of the parameter 
space. If R is not big and mutation is high enough, the 
lower classes contribute mostly to higher ones, get less 
populated and, consequently, are more sensible to fluc- 
tuations. This scenery seems appropriate for an abrupt 
loss of lower classes and that is exactly what results from 
an empirical study of some combinations of parameters, 
that the ocurrence of the anomalous behavior in Muller's 
ratchet is favored by high mutation rates and low fertil- 
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FIG. 1: Dynamical behavior of the average mutational load 
of the least loaded class for R — 2 and U = 0.6. At t = 15, 
from top to bottom, s = 0.1, 0.05, 0.03, 0.0, 0.2 and from 
this value s increases to 0.7 in steps of 0.1. The anomalous 
activity of Muller's ratchet found by FCH appears clearly. 
The continuous lines just aid visualization of the discrete time 
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FIG. 2: Dynamical behavior of the average mutational load 
of the least loaded class for R — 2 and U = 0.2. At i = 15, 
from top to bottom, s = 0.0, 0.03, 0.05, 0.1 and from this 
value s increases to 0.7 in steps of 0.1. Here, Muller's ratchet 
presents a monotonic dependence on s. The lines are as in 

pig. m 



ity For instance, Fig. [21 where the mutation rate is just 
a bit lower than in Fig. [TJ already shows a monotonic 
response of Muller's ratchet to variation in s. 

Finally, the SBP theory gives the condition for the 
halting of Muller's ratchet also under decay. The rea- 
soning is analogous to the neutral case. By definition, if 
the ratchet does not halt, the minimum number of mu- 
tations in the system grows unrestrictly. In this case, 
since (1 — s) J decreases monotonically to zero with the 
mutational load, no matter how small is s > 0, there is 
a finite time when the minimal mutational load j* is so 
high that the average number of perfect replicates of each 
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of the fittest individuals, Re~ u (l — s) J , is smaller than i 
and even the fittest subpopulation certainly goes extinct, 
and obviously extinction is the fate of the whole popu- 
lation. Thus, no population can stand endless mutation 
accumulation, and the ratchet is certainly halted in any 
indefinitely surviving lineage. Besides that, Re~ u > 1 
turns out to be the necessary condition for survival un- 
der decay, since the master class is unaffected by s and 
performs better than all other classes in surviving. 

Hence, whatever be s, the halting criterium of Muller's 
ratchet is Re~ u > 1, although it does not assure the sur- 
vival of the population, ft is noteworthy this result is a 
direct consequence of the basic theory of branching pro- 
cesses. But this formalism also gives elaborate quantita- 
tive tools that allow a thorough study of models of mu- 
tation accumulation in growing lineages. In particular, 
any decay law can be analyzed and the assumptions con- 
cerning the asymptotic survival of individuals with few 
mutations now can be evaluated at low computational 
cost. 

Besides that, the stochastic dynamics of a growing lin- 



eage described in the present paper, as well as analogous 
dynamical solutions recently found for the deterministic 
behavior of infinite po pula tions evolving on multiplica- 
tive [HI and truncated |l9l] fitness landscapes, may prove 
useful in the construction of general theoretical models 
(an early example is [8() suited for describing popula- 
tions mixing growth with bottlenecks, like that evolved 
under the serial transfer protocol of experimental evolu- 
tion 0, m 

These important themes are out of the scope of this 
short paper and will be discussed elsewhere. 
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